Can Cosmic Structure form without Dark Matter? 
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One of the prime pieces of evidence for dark matter is the observation of large overdense regions in 
the universe. Since we know from the cosmic microwave background that the regions that contained 
the most baryons when the universe was ~ 400, 000 years old were overdense by only one part in 
ten thousand, perturbations had to have grown since then by a factor greater than (1 + z„) ~ 1180 
where z* is the epoch of recombination. This enhanced growth does not happen in general relativity, 
so dark matter is needed in the standard theory. We show here that enhanced growth can occur 
in alternatives to general relativity, in particular in Bekenstein's relativistic version of Modified 
Newtonian Dynamics (MOND). The vector field introduced in that theory for a completely different 
reason plays a key role in generating the instability that produces large cosmic structures today. 



Introduction. Dark Matter was introduced long ago to 
explain galactic rotation curves [jj and large velocities 
in galaxy clusters 0). Over the past decade the case for 
dark matter has gotten stronger: galactic rotation curves 
still appear to diverge significantly from what is expected 
from observations of visible matter; gravitational lensing 
is able to map the mass distribution in a galaxy or cluster 
and this mass distribution often does not coincide with 
the luminous matter and large overdense regions in 
the universe can be explained only if dark matter was 
around early on to seed structure formation when the uni- 
verse was of order several hundred thousand years old. At 
the same time, the case that dark matter does not con- 
sist simply of protons and neutrons which do not emit 
light has also gotten stronger 0, ■ To explain the as- 
tronomical and cosmological observations, therefore, we 
apparently need to introduce a new fundamental particle 
which has not yet been observed in accelerators. 

There is one way of avoiding this conclusion: perhaps 
the implicit assumption that gravity is described by gen- 
eral relativity is incorrect. Perhaps a fundamental theory 
of gravity which differs from general relativity on large 
scales can explain the observations without recourse to 
new, unobserved particles. Now more than ever before, 
there are very good reasons to explore this idea of modi- 
fying gravity. For, the case for dark energy also hinges on 
the assumption that general relativity describes gravity 
on large scales. Dark energy is even more difficult to ex- 
plain in the context of fundamental theories than is dark 
matter, so it seems almost natural to look at gravity as 
the culprit in both cases. 

Perhaps the most direct piece of evidence for dark mat- 
ter is from galactic rotation curves, which are much flat- 
ter than the i? _1//2 fall-off expected from observations 
of visible matter. Over 25 years ago, Milgrom pro- 
posed a MOdified Newtonian Dynamics (MOND), which 
diverges from Newtonian theory when the gravitational 
acceleration is less than ao ~ (200krnsec _1 ) 2 /(10kpc). 



Argument has raged for two decades as to whether this 
modification is consistent with a wide variety of observa- 
tions 0, 0, Hol HH n~3 | . Part of the difficulty in assessing 
MOND is that it makes no claims to be a comprehensive 
theory of gravity, so, for example, it is impossible to test 
it with cosmological observations. 

Recently, Bekenstein constructed a fully relativis- 
tic covariant theory which reduces to MOND in the ap- 
propriate (static) limit. In addition to the gravitational 
metric (a tensor field), Bekenstein's theory contains a 
scalar field and a vector field; hence he called it Ten- 
sorVectorScalar, or TeVeS. We can now compare the pre- 
dictions of TeVeS to those of general relativity. There are 
two features of the new theory that are particularly inter- 
esting and seem worthy of further study. First, several 
new constants must be introduced into the TeVeS La- 
grangian and one of these is of order ao, the MOND scale. 
This is no surprise of course since TeVeS is designed to 
reduce to MOND in certain limits. The interesting part 
is that ao is roughly the same order of magnitude as the 
fundamental scale introduced in quintessence or modi- 
fied gravity theories designed to explain the acceleration 
of the universe. The motivation for TeVeS had nothing 
to do with the cosmic acceleration. Is it just a coinci- 
dence then that the fundamental scale needed is the one 
required to obtain acceleration? Or is this an indication 
that we are on the right track in our quest to explain 
away the dark sector by modifying gravity? 

The second intriguing aspect of TeVeS, and the one 
we focus on h ere, was recently uncovered by Skordis and 
collaborators |l5j . To understand the importance of this 
feature, it is necessary to state an obvious hurdle that 
any no-dark-matter theory must overcome. This hurdle 
is best depicted in Figure which shows the observed 
power spectrum of matter in the universe in a theory 
without non-baryonic dark matter. This theory fails to 
describe observations in two ways: (i) the shape of the 
power spectrum is off and (ii) the amplitude is far too 
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FIG. 1: Power spectrum of matter fluctuations in a the- 
ory without dark matter as compared to observations of the 
galaxy power spectrum. The observed spectrum [T^l does 
not have the pronounced wiggles predicted by a baryon-only 
model, but it also has significantly higher power than does 
the model. In fact A 2 , which is a dimensionless measure of 
the clumping, never rises above one in a baryon-only model, 
so we would not expect to see any large structures (clusters, 
galaxies, people, etc.) in the universe in such a model. 



small. The first failure has been exploited by many au- 
thors to prove the existence of non-baryonic dark mat- 
ter [lH ITtI , the statistical significance for which now 
exceeds 5-sigma. The second failure is often ignored be- 
cause analysts typically marginalize over the amplitude 
of the power spectrum on the grounds that the power 
spectrum of galaxies is likely to differ by an overall nor- 
malization factor (the bias) from the power spectrum of 
matter. But a baryon-only model fails miserably at get- 
ting anywhere near the amplitude required to generate 
galaxies and galaxy clusters even with an absurd amount 
of bias. So if we really want to do away with dark matter, 
we need to find a mechanism of growing perturbations 
faster than in standard general relativity. This is pre- 
cisely what Skordis et al. 0, 0] seemed to have found 
in their treatment of perturbations around a smooth cos- 
mological solution in TeVeS. Here we aim to move beyond 
their numerical treatment to isolate what is causing en- 
hanced growth. Our motivation goes beyond TeVeS, as 
the exact Lagrangian in will almost certainly need to 
be altered even if the general idea turns out to be correct. 
Indeed, as shown in Fig^ even if structure grows faster 
than in the standard theory, the shape of the baryon- 
only spectrum does not match the observations. Rather, 
we want to understand generally how to modify gravity 
such that it solves not only the galactic rotation curve 
problem but also the cosmological structure problem. 
Cosmology in TeVeS. Ordinary matter couples to the 
gravitational metric <? M „ in the standard way in the TeVeS 
model. The metric which couples to matter, though, does 
not appear in the standard way in the Einstein-Hilbert 
action. Rather, it is useful to define a new tensor g^ v 



which is a functional of g^ v and a scalar field <f> and a 
vector field A^. Specifically, 
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defines g^. The action of <? M „ is the standard Einstein- 
Hilbert action. The scalar and vector fields have dynam- 
ics given, respectively, by the actions S s and S v : 
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where /i is an additional non-dynamical scalar field, 
F^ u = A^ v — A VtfJj , and indices are raised and lowered 
with the metric g M „. The potential V(fi) is chosen to 
give the correct non-relativistic MONDian limit. We will 
consider the form proposed by Bekenstein [13j: 
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[A(4 + 2/t - 4/x + /i 3 ) + 2 In (fi - l) 2 ] 



(3) 

with fi = fi/fiQ. There are three free parameters that 
appear in the TeVeS action: (1q, £b and Kb- The pa- 
rameter A in the vector field action is completely fixed 
by variation of the action. 

Armed with this action, we can solve [l3l Hif for 
the evolution of the scale factor a of a homogeneous 
Friedman-Robertson- Walker (FRW) metric. This evo- 
lution turns out to be very similar to the standard case, 
with several small deviations. First, Newton's constant 
gets generalized to Ge~ 4( ^/(1 + d(f>/d\n(a)) 2 . Second, the 
Friedman equation governing the evolution of a has, in 
addition to the standard source terms of the matter and 
radiation energy densities, the energy density of <f>: 
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FIG. 2: Evolution of homogeneous TeVeS fields. Dashed line 
shows logarithmic approximation for cf> valid in the regime 
when n is constant. In that regime, scales as the ambient 
density, with the ratio equal to (6po)~ in the matter era. 
Early on, p^/ptotai = -</>= 15/(4/i). 

The TeVeS modifications to the standard cosmology 
then depend on the evolution of the scalar field 4>. Dur- 
ing the radiation dominated era, is much smaller than 
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the dominant radiation density, with the ratio growing as 
a 4 / 5 0|- As shown in Fig.|3 once a gets large enough, 
fi becomes constant and <f> ~ constant + ln(a)/2p Q . The 
energy density in cj) then scales as a -3 exactly like mat- 
ter, with the ratio fixed at (6/tto) . For a variety of 
reasons ^jj, large values of po(~ 100 — 1000) are pre- 
ferred, so pif, is much smaller than the other densities at 
all times and the effective Newton constant differs at only 
the percent level from the standard one. Homogeneous 
expansion is there the same in TeVeS as in standard gen- 
eral relativity. 

Perturbations in TeVeS. We start the perturbation ex- 
pansion by writing the metric which couples to matter 
as 

9oo(x,t) = -a 2 (r)(l-2*(f,r)) (5) 
5y (f,r) = o 2 {t){\ + 2^{x,t))5 ij (6) 

where r is conformal time; similarly we expand the den- 
sity of the matter (and/or radiation) field as 



p(x,r) = p(r) (1 + S(x,t)) . 



(7) 



To consider the evolution of perturbations in TeVeS, 
we need to perturb not only the matter and the FRW 
metric fields, but also the new fields: 
The scalar perturbation can be written as 



and mint. 



(x,t) = <f>(r) + ip(x,r) 



(8) 



where <f> is the zero order field introduced above. 

The vector field requires a little more thought • 
In general, perturbations to a vector field A^ will be 
described by four independent functions. In this case, 
though, the two transverse spatial components decou- 
ple from the set of scalar perturbation equations, so we 
can neglect them. Further simplification follows from the 
fact that the vector field in TeVeS is subject to the con- 
straint M'A^ ee g^A^A,, = -1. This fixes the time 
component of the perturbation, so we need track only 
the longitudinal component of A^. In detail, the zero 
order A^ can be chosen to have only a time component. 
The constraint then sets that time component to ae~^ 
(since g 00 = —a~ 2 e 2 ^), so the perturbed vector field can 
be written 



where A^ = (1, 0, 0, 0) and 



= (\& + tp, a) 



(9) 



(10) 



That is, the time component of the perturbation is con- 
strained to be a combination of the perturbations to the 
metric ('&) and the scalar field (ip). Skorids et al. 
called the longitudinal component of the perturbation 
a (with no index); specifically, Va ee a or equivalently 
V 2 a ee V • a. 



The perturbations to the metric, matter, radiation, 
and TeVeS fields are governed by a set of coupled differ- 
ential equations. Ref. suggested that perturbations 
in the scalar field may induce enhanced growth in the 
matter perturbations. We have found [l9| that this is 
not so: the scalar perturbations are small early on and 
then oscillate about Q/p,o in the matter epoch. This is 
far too small a value to impact density perturbations. 

Rather we find that the vector perturbations are 
the key to enhanced growth. The equation governing 
vector perturbations is: 



a + b 1 a + b 2 a = S[$, 



(11) 



where the source term on the right is a functional of the 
metric perturbations and the coefficients on the left are 
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Here ai is defined as 
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There is one limit in which analytic solutions exist for 
the homogeneous part of Eq. (|llf> : when the background 
quantities are such that <j) and <fi vanish. This is a fairly 
good approximation because we know that <f> has little 
impact on the zero order dynamics. In this limit, b\ — > 
2a/ a and bi — > a/ a. For simplicity then consider a matter 
dominated universe so that b\ — 4/r and 62 = 2/r 2 . For 
reasons that will become clear soon, let us write 
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In the limit we are now considering, e = 0. 

The homogeneous solutions to Eq. i|ll|) scale as r p with 
the powers determined by solving an algebraic equation 
so that 
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P± = — ±2^9-8(1 + 4 



(15) 



So when e = 0, the two homogeneous modes scale as r~ 2 
and r _1 . The particular solution will then dominate; in 
the matter era, the dominant source term is — so 
a = — The top panel of Fig. [3] shows that, when 

Kb is not too small, a does indeed follow the particular 
solution. 

However, when Kb is small, the term in Eq. H12|) mul- 
tiplied by Kg 1 cannot be neglected. The ratio of this 
term to d/a is defined as e. Both terms in square brack- 
ets scale the same way, but the first term is quite a bit 
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larger. In a matter dominated universe, 8irG(p + P)a 2 
3a 2 H 2 — > 12/r 2 , so using the logarithmic solution for . 
we find 



-12 In (a/5 x 1(T 5 ) 



(16) 



So as long as Kb^o is of order 100 or smaller, the coeffi- 
cient of a in Eq. (|11|) . 1 + e, will eventually become neg- 
ative, signifying enhanced growth. Note that because of 
the slow, logarithmic growth of 0, the enhanced growth 
will not be exponential. It will, however, be a grow- 
ing mode, which deviates from the particular solution 
— \I/t/3. The top panel of Figure 01 shows that the qual- 
itative aspects of this analytic solution do emerge in the 
full numerical results. 
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FIG. 3: Evolution of cosmological perturbations (unnormal- 
ized) in a TeVeS model with Q m = 0.3 (baryons only), 
Ho — 200 and low value Kb = 0.07 (dashed) and a high 
value of Kb = 1 (dot-dashed). Top panel shows that vec- 
tor perturbations become unstable for the low value of Kb- 
Solid curve is the particular solution a = — tyr/3. Second 
panel shows that this induces a large difference between the 
two Newtonian potentials. Third panel shows that this drives 
enhanced growth in the density perturbations as compared to 
standard ACDM (solid curve) if Kb is small; density pertur- 
bations in the large Kb case are smaller than in ACDM due 
to the absence of dark matter. In all cases the wavenumber 
is k = 0.5Mpc _1 

The growing vector field drives the two Newtonian po- 
tentials to differ from one another as seen in the middle 
panel of Fig|3] Recall that in general relativity, this dif- 
ference is sourced only by anisotropic stress. In TeVeS, 
the vector field also sources the difference 01- This dif- 
ference in turn drives enhanced growth in the density per- 
turbations as shown in the bottom panel, precisely the 
kind of growth needed to generate large structures from 
the small inhomogeneities present at recombination. We 



have verified that this growth does not occur if vector 
perturbations are turned off. 

How generic are these ideas of vector instabilities and 
their subsequent impact on the two Newtonian poten- 
tials? Several authors [2JJ, have pointed out that, 
in the context of the general Lagrangian studied in 
[22| . vector instabilities exist for a wide range of coeffi- 
cients. So vector instability generally seems quite plausi- 
ble. Bertschinger [23| has pointed out that the key input 
from modified gravity models is the source term for the 
difference between the two Newtonian potentials. So it 
is not surprising that this difference plays an important 
role in TeVeS. It is possible that the enhanced growth dis- 
cussed here will emerge naturally from a modified gravity 
model which explains the acceleration of the universe. 

Conclusions. Motivated by observations of gravitational 
lensing, Sanders Q was the first to introduce a vec- 
tor field into the MOND framework. Bekenstein's re- 
cent theory 0] elevated this field to be dynamical. We 
have shown here that, at least in Bekenstein's framework, 
the vector field can also source enhanced growth in the 
cosmic density perturbations. This is by no means the 
last word on confronting cosmological data with MON- 
Dian theories, but it does overcome perhaps the primary 
cosmological hurdle faced by a baryon-only model. The 
enhanced growth enables the very small perturbations 
at recombination to grow into the large structures we 
see today. Among the significant problems that remain 
are the need to match the observed galaxy power spec- 
trum on large scales and the well-measured series of peaks 
and troughs in the CMB spectrum and the apparent mis- 
match between mass and light in galaxy clusters 0. 
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